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Abstract. In this note we improve an algorithm from a recent paper 
by Bauer and Bennett for computing a function of Erdos that measures the 
minimal gap size f(k) in the sequence of integers at least one of whose prime 
factors exceeds k. This allows us to compute values of f(k) for larger k and 
obtain new values of f(k). 
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1 Introduction 



For any integer m we let P(m) be the largest prime factor of m with the 
convention P(0) = P(±l) = 1. Let H n) k be the product of k consecutive 
integers, starting with n, i.e. 

n n ,fc = n(n + l)--(n + k — l). 

By a theorem of Sylvester (see [U]) Yi n k is divisible by a prime p > k 
whenever n > k, and thus, following Erdos [7], we define f(k) to be the least 
integer with the property that 

p (n nj/(fc) ) > k, 

if n > k. Standard heuristics for the size of gaps between consecutive primes 
lead one to expect that order of magnitude of f(k) is (log A;) 2 . 
The following table gives known values of f(k): 



k f(k) 


k f(k) 


k f(k) 


1 1 

2 2 
3-4 3 
5-12 4 


13-40 6 
41-46 7 
47- 58 8 
59 - 60 9 


61-113 14 
114 13 
115-150 12 
151 - 178 14 



The values of f(k) for k < 10 were computed by Utz [TB] and extended to 
k < 42 by Lehmer [llj, to k < 46 by Ecklund and Eggleton [4], to k < 73 
by Ecklund, Eggleton and Selfridge [5], [6] and finally to k < 178 by Bauer 
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and Bennett |1J. Bauer and Bennett in the same paper also disproved an 
assertion of Utz that / is monotone. 

In this paper we compute the values of f(k) for k < 268. Our results can 
be summarized in the following theorem. 

Theorem 1. For 179 < k < 268 the values of f(k) are as follows: 





k 


/(*) 


179 


- 222 


14 


223 


-268 


16 



Note that Bauer and Bennett, although greatly extending the set of k 
such that f{k) is known, by 105 values of k, did not find any values of f(k) 
such that f(k) > 14. In Theorem 1 we find the new largest proven value of 
f(k), the first after nearly 40 years (the previous being [6]). 

2 The algorithm 

Lemher [10J searched for two consecutive smooth integers, satisfying P(z(z + 
1)) < Pt, where p t is the t-th prime. One can write x = 2z + 1 and see that 
P(z(z + 1)) < p t iff P(x 2 — 1) < p t . Writing x 2 — 1 = dy 2 , where d is 
squarefree, leads to the Pell equation 

x 2 - dy 2 = 1, (1) 

where P(dy) < p t . This implies that the solutions x + y\/d we are searching 
for are powers of the fundamental solution u + v \/d, i.e 

x + yVd = (it + v\/d) n . 

By classic results on primitive divisors (see [3]) and the work of Lehmer (see 
[9]), if we want P(y) < pt, then 

. f pt + l 10 \ 
n < max < , 12 > . 

Thus one needs to consider only finitely many Pell equations and for each 
one only some of the first solutions. 

After finding all pairs of smooth consecutive integers in this way, one can 
search through the results and find larger strings of consecutive integers. 

Bauer and Bennett [1J improve on this strategy by the following clever 
argument: for fixed m < 3 and t suppose we are searching for all integers n 
satisfying 

P(n n , m ) < p t . (2) 
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We can split the indices 0, . . . , m — 1 into |_xJ + disjoint pairs 

(i,i + 2), where i = 0, 1 (mod 4). Set t = 7c(m — 1). By the Dirichlet 
principle, we can find an index % such that in + i)(n + i + 2) is divisible by 
at most 

t-t 



N 

of the primes from the set 



LfJ + L^J 



{Pt +i,- ■ ■ ,Pt}- (3) 
Now one writes X = n + % + 1 and Y = (n + z)(n + z + 2), and gets 



x 2 - L>y 2 = l, 



where D is squarefree and divisible only by some of the first to primes and at 
most N of the primes from the set Then one proceeds exactly as Lehmer, 
with the difference that in this approach P(X) < pt also has to hold. This 
lowers the number of Pell equations one needs to consider from 2* — 1 to 



and equally important, reduces the size of the Pell equations. 

The bottleneck of both the algortihms of Lehmer and Bauer and Bennett 
is solving the Pell equation. The main difficulty in solving the Pell equation 
is the size of the solutions, as it grows exponentially with respect to the size 
of the coefficient d form ([T]). This means that just writing down the solution 
in standard representation takes exponential time. In [12J, Luca and the 
author used compact representations of the solutions to the Pell equations to 
overcome this difficulty and managed to extend Lehmer's results (see also [H] 
for another application of this approach). Note that the algorithm described 
in [T2] is still the best up to date if one wants to find 2 or 3 consecutive 
smooth integers. 

A compact representation of an algebraic number f3 £ Q(Vd) is a repre- 
sentation of of the form 

k 

where ctj = [aj-\-bj\fd)/cj, where dj,bj,Cj,k G Z, and there exist polynomial 
(in log d) upper bounds for k, ctj, bj and Cj. A detailed description of compact 
representations and their use can be found in [8]. Using compact represen- 
tations, the Pell equation is solved in two steps. First, the regulator of the 
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appropriate quadratic field is computed. Once the regulator is known, one 
computes a compact representation, in polynomial time. By far the harder 
part of this process is the computation of the regulator. The fastest uncon- 
ditional algorithms are still exponential with respect to d, while the fastest 
known algorithm, Buchmann's subexponential algorithm (see j2]), depends 
on the Generalized Riemann Hypothesis. 

For our purposes, the only algorithm fast enough is Buchmann's algo- 
rithm. We will perform a check, using continued fractions, that will for each 
case tell us either that the output of Buchmann's algorithm is unconditionally 
correct or that we can disregard this case. This implies that, although we 
cannot unconditionally solve the Pell equation in subexponential time, we 
can determine unconditionally whether it has a smooth solution (and find 
any smooth solutions) in subexponential time. 

We do this by the following method, first used in [12]: suppose we want 
to find the solutions of x 2 — dy 2 = 1, where x and y are n-smooth (d is a 
fixed n-smooth integer). We first compute the regulator R d of Q(y/d) using 
Buchmann's subexponential algorithm. As mentioned, the correctness of the 
output Rd is dependent on the GRH, but what is unconditionally true is 
that the output is equal to mRa, where m is some positive integer. From the 
obtained value we construct, in polynomial time, a compact representation of 
the m-th power of the fundamental unit, and from it a compact representation 
of the m-th solution x m + y m ^fd of our Pell equation. We define z to be the 
n-smooth part of y m , which is computed using the algorithm for performing 
modular arithmetic on compact representations from [T3] . 

Let x\ + yiVd be the fundamental solution of our Pell equation and 
Xk + ykVd = (x\ + y\\fd) k . If the GRH is false and Buchmann's algorithm 
does not return the exact value of the regulator this will mean that it might 
happen that we will miss some smooth solutions of our equation. Suppose 
this happens. This means that m > 1 and yk is n-smooth for some k that 
is not a multiple of m. But we have yi\yk, so y\ is also n-smooth, implying 
y\\z (because yi\y m )- This obviously implies y\ < z. So if something has 
gone wrong, we will find the solution yi that is smaller or equal to z and 
strictly smaller than y m . We can use continued fractions to check whether 
such a yi exists. If we do not find such a y\ (and indeed we never do, as it's 
existence would disprove the GRH), it means either m = 1 or m > 1 but we 
are not missing a smooth solution. In any case, this removes the dependence 
of our results on the Generalized Riemann Hypothesis. Note that although 
y m does become very large, z in practice does not. This means that it is 
usually, for our purposes, enough to compute the first 10 or so convergents of 
the expansion of y/d to a continued fraction, to prove that the hypothetical 
yi does not exist. 



4 



Using compact representation, in combination with this check, cuts down 
the space needed from exponential to polynomial with respect to d and the 
time needed to solve the Pell equation from exponential to subexponential. 

In essence, in this paper we combine the clever approach of [1] and the 
powerful methods for finding smooth solutions of the Pell equations from 

P31- 

It is worth noting the somewhat counter-intuitive implications of the 
mentioned algorithms: it is easier to solve the Pell equation (using compact 
representations) than it is to write down the solution in decimal notation, and 
it is easier to unconditionally find all smooth solutions of the Pell equation 
or determine that there are none (this takes subexponential time) than it is 
to unconditionally solve the Pell equation (this takes exponential time). 

3 Results 

We run our algorithm for two pairs of parameters (m,t), these pairs being 
(47, 14), and (56, 16). Note thet p 47 = 211, p 56 = 263 and P(n 318>13 ) = 163, 
■P(ni330,i5) = 223, implying 

f(k) > 14 for 178 < k < 223 and f(k) > 16 for 224 < k < 268 (5) 

By proving that there is no integer n > k for 178 < k < 222 satisfying 
-P(II rai i4) < k implies f(k) = 14 for 178 < k < 222. In the same way one 
proves that f(k) = 16 for 223 < k < 268. 
The case (m, t) = (47, 14) makes us solve 




= 342948991 



Pell equations, while the case {m,t) = (56, 16) makes us solve 
M 2 = -1 + 2 6 ( ■ ) = H68680703 

Pell equations. 

We obtain that there do not exist such integers n, thus proving Theorem 1. 
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